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ABSTRACT 

The distributions of deuterated molecules in protoplanetary disks are ex- 
pected to depend on the molecular formation pathways. We use observations 
of spatially resolved DCN emission from the disk around TW Hya, acquired 
during ALMA Science verification with a ~ 3" synthesized beam, together with 
comparable DCO+ observations from the Submillimeter Array, to investigate dif- 
ferences in the radial distributions of these species and hence differences in their 
formation chemistry. In contrast to DCO+, which shows an increasing column 
density with radius, DCN is better fit by a model that is centrally peaked. We 
infer that DCN forms at a smaller radii and thus at higher temperatures than 
DCO"*". This is consistent with chemical network model predictions of DCO"*" 
formation from H2D"'" at T<30 K and DCN formation from additional pathways 
involving CH2D"'" at higher temperatures. We estimate a DCN/HCN abundance 
ratio of ~ 0.017, similar to the DCO"'"/HCO'^ abundance ratio. Deuterium frac- 
tionation appears to be efficient at a range of temperatures in this protoplanetary 
disk. These results suggest caution in interpreting the range of deuterium frac- 
tions observed in Solar System bodies, as multiple formation pathways should be 
taken into account. 

Subject headings: astrochemistry — circumstellar matter — planetary systems 
— molecular processes — planet-disk interactions 
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Introduction 



Deuterium fractionation in interstellar cloud cores, protostars and Solar System bod- 
i es is frequeritly us e d to infer important aspects of their physical and chemical histories 



( iBrown et al.l Il988l : iHerbst fc van Dishoeck 



20091 : iMumma fc CharnlevI 1201 If ). The deu- 



terium enhancement in the Earth's sea water, with respect to the cosmic abundance, has 
been used as an argument for volatile delivery to the young Earth from cold, comet-like 
planetesimals, even though most presen t-day comets have about a factor two higher D/H 
ratio compared to the Earth's sea water (IMumma fc CharnlevI 1201 111. Recent observ ations of 
Earth-like deuterium fractionation in a Jupiter family comet ( iHartogh et al.l 120111 ) has lent 
further support to this scenario. Deuterated molecules, notably DCO^ and DCN, are also 
detected in a handful of protoplanetary disks, the sites of ongoing pla net formation, with 



aburidances that imply deuterium enhancern e nts by orders of magnitud e (Ivan Dishoeck et al. 



2003 



Guilloteau et al. 


2006; 


Oi et al. 


2008; 


Obere; et al. 


2010. 


2011) 



Understanding how 

and under what conditions these deuterated molecules form in disks is key to use deuterium 
fractionation to probe physical conditions, and to explain the deuterium fractionation pat- 
tern in the Solar System. The formation chemistry of DCN is of particu lar interest since the 
DCN/HCN ratio can be measured in both comets (e.g. iMeier et al.lll998l ) and disks, enabling 
direct comparison. 

Deuterium levels in molecules are enhanced above the cosmic deuterium abundance of 
~ 2 X 10~^ because of differences in zero-point energy between deuterated species and the non- 
deuterated equivalents. The reverse reactions are efficiently inhibited at low temperatures, 
resulting in deuterium fractionation. However, the "low temperature" regime is different for 
different reactants. Deuterated species that form from reactions wit h H^D"*", in p articular 
DCO"*", are expected to be abundant only at temperatures below 30 K (IWillacvll2007l) because 



of th e modest activation energy of 230 K for converting H2D+ back into ( ICerlich et al. 



2OO2I ). The activation ener gy barrier for the conversion between CH2D"'" and CHjj' is higher, 
390 K (lAsvany et al.ll2004l ). and CH2D+ can therefore drive a warmer deuterium chemistry. 
Theoretically, DCN may form through reactions with the CH2D+ reaction product CHD, as 
well as with H2D+, DC O"*", D, and through grain surface rea ctions followed by thermal or 
non-thermal desorption (lAikawa &: Herbstll200ll ; IWillacyl 120071 ). Observations of warm DCN 
in the Orion bar provides observational evide nce for that the C H2D"'" reaction pathway is 
important under some interstellar conditions (jParise et al.ll2009l ). The formation pathways 
of DCN that are efficient in disks have yet to be constrained by observation. 

TW Hydrae is the most well-studied protoplanetary disk because of its proximity (51 pc) 
and near face-on viewing geometry that allows for direct investigation of the radial distri- 
bution of molecular emission. Millimeter emission from CO, HCO+ 3-2, DCO"*" 3-2, HCN 
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3-2 and DCN 3-2 lines have been spatially resolved using the Submillimeter Array (SMA). 
Analysis of DCO+ emission showed that the DCO"*" column density increases with radial dis- 
tance from t he star, out to ~ 90 AU, consistent with the low-temperature formation pathway 
from H2D"'" (iQi et al.ll2008[ ). The DCN emission was too weak for detailed modeling, but 
newly released ALMA science verification data on DCN 3-2 towards TW Hydrae have higher 
signal-to-noise and allow for significant constraints. In this paper, we use these new data 
to compare the DCN and DCO+ distributions, and to test if these data suggest different 
pathways to deuterium enhancements in these species as predicted by some models of disk 
chemistry. 



2. Spatially Resolved Observations of DCN and DCO+ towards TW Hya 

TW Hya was observed on April 20th and 23rd, 2011 in ALMA band 6 as a part of ALMA 
Science Verification. The calibrated visibilities and CLEANed reference images were released 
and made available to the public in August 20110. The.e oteervations include the DCN J=3- 
2 line at 217.238 GHz, which was acquired with a channel spacing of 0.16 km s~^ (resampled 
to 0.2 km s"^) and synthesized beam size 2'.'8 x 2'.'3, with an rms of 12 mJy beam~^ in the 
line-free chann els. These data are comparable in resolution to the previous observations of 



Qi et al.l ( 120081 ). with a noise level 14 times better for the ALMA data compared to the SMA 



observations. 

Figure [1] displays the integrated DCN 3-2 emission toward TW Hya. The DCN emission 
is peaked on the TW Hya stellar position and spatially compact, barely resolved with the 
ALMA beam. For comparison. Figure [T] shows the integrated DCO^ 3-2 emission, which 
does not peak at the stellar position and appears more extended. Using the DCO+ velocity 
channel maps, especially the peanut-shaped central one that provides strict constraints on 



the DCO+ radial distribution, iQi et al.l ( 120081 ) demonstrated that the DCO"*" column density 



increases with distance from the central star out to 90 AU. Since the excitation characteristics 
of the DCN and DCO+ J=3-2 lines are similar, the differences in the Fig. 1 images suggest 
different abundance distributions. In particular, DCN appears to be more abundant in the 
inner (warmer) disk regions compared to DCO+. 



ht tps : / / almascience .nrao.edu/ almadat a/ sciver /T WHyaBand6 / 
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3. The DCN radial profile 



T o constra i n the radial distribution of the DCN column density, we use the same meth- 
ods as lQi et al.l ( 120081 ) to model observations of DCO"'", HCO+ and HCN emission from TW 
Hya. First, we ad opt a physical structu re of the TW Hya disk with density and temperature 
distributions from lQi et al.l (120041 12006| ). Second, we assum e DCN ex i sts in a disk layer with 
vertical boundaries given by the best-fit values for HCN in iQi et al.l (120081 ): the boundaries 
of DCN are not very well determined from the present data a nd chemical models predict a 
similar vertical distribution of DCN and HCN (jWillacyl l2007l ). Third, we model the DCN 
column density in this layer as a power law, A^io x (r/10)", where A^io is the column density 
at 10 AU, r the distance from the star in AU, and a the power-law index. Based on the DCN 
data, the model is cut off at 100 AU, which is in agreement with previous models of DCO"*" 

and HCN. For HCN and DCO+, lOi et all toO^ found a 1 and a ~ 2, respectively. 

We therefore tested models of DCN column density with three different values of power-law 
index, a = —1, 0, and 2, to investigate if the DCN distribution can be modeled similarly 
to DCO^, or if it can be better described by a fiat or decreasing function of radius. A 
more detailed characterization of the shape of the DCN radial profile requires higher spatial 
resolution than currently available. 

For each power-law model, we optimize A'^io by minimizing the value, the weighted 
differences betwee n the observed a nd modeled complex visibilities. The techniques are de- 
scribed in detail in lQi et al.l (120081 ) . and the best-fit models shown in Figure [2] and Figure [31 
Unlike the DCO"*", the DCN J=3-2 line is complicated by hyperfine structure. The cen- 
tral DCN line quartet at 217.23863 GHz is separated by 0.23 MHz (0.32 km s"^) from a 
singlet at 217.23840 GH:J^, with 22% of the intensity of the quartet. This is close enough 
in frequency and strong enough in relative intensity to affect the shape and strength of 
emission in the channel maps. We use the Monte Carlo radiative transfer code RATRAN 



(IHogerheijde &: van der Takll2000l ). taking into account the DCN hyperfine structure, to pro- 
duce model visibilities that are sampled at the appropriate spatial frequencies for comparison 
with the ALMA data. The relative population of the hyperfine levels are assumed to be in 
LTE, and the relative optical depth of the hyperfine transitions including line overlap is 
accounted for accurately 

Figure [3] shows a comparison between the observed DCN J=3-2 line channel maps and 
the best fit models for power-law indices a = —1, and 2. The model with a = —1 clearly 
provides the best fit to the data, as models with a = and a = 2 underestimate the emission 
in the line wings, show offsets from the observed peaks, and appear more elongated in the 



^http://spec.jpl. nasa.gov/ 
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central channel than the data. Note that the peanut-shaped f eature in the sy stemic velocity 
channel map of DCN Model 3 is not as obvious as for DCO"^ ( IQi et al.ll2008[ ) because of the 
blending of DCN hyperfine components. The better fit of the model with a negative power 
index is reflected in the values of the models, which increase from 2296194 for Model 1, 
to 2298334 for Model 2, and to 2300603 for Model 3. (Note that we have not excluded the 
possibility that DCN is even more centrally peaked than assumed in Model 1.) 



Using Model 1 and the HCN profile reported by iQi et al.l (120081 ). we calculate the 
DCN/HCN abundance ratio to be ~ 0.017 in the TW Hya disk. This result does not de- 
pend strongly on the model power-law indices; varying the power-law indices by ±2 changes 
the ratio by less than 30%. This ratio is simil ar to the previously meas ured DCO"'"/HCO^ 
abundance ratio of 0.035 in the TW Hya disk (Ivan Dishoeck et al.ll2003l ). 



4. Discussion 

The observations demonstrate that DCN is centrally peaked on the size scales observed, 
ruling out common radial distributions for DCN and DCO"*". The presence of DCN closer to 
the star, in regions that are warmer due to stellar irradiation, together with the comparable 
estimates for the DCN/HCN and DCO^/HCO^ abundance ratios, suggests a mechanism 
for efficient deuterium enhancement in disk material at higher temperatures than implied 
by reactions solely with H2D+ at T<30 K. This agrees with standard model predictions, in 
that DCN can form from gas phase reactions involving CH2D"*" at T>30 K. However, there 
are several alternative chemical mechanisms that have the potential to affect the relative 
distributions of DCO"'' and DCN, including photochemistry, additional gas-phase reactions 
that result in HCN deuterium fractionation, grain surface formation of DCN followed by 
desorption, and molecule specific destruction such as depletion onto interstellar grains. Since 
the origins of DCO"*" are reasonably well understood, we will focus on the mechanisms that 
affect the DCN distribution. 

Before doing so, we note that it is possible for dynamical processes to transport cold 
chemistry products to warmer disk regions. If the destruction rates are slow compared to 
formation, then high deuterium fractionation at T > 30 K does not a priori imply warm 
deuterium chemistry. In the context of the DCN and DCO+ distributions observed in the 
TW Hya disk, both DCN and DCO"*" could form cold, and then DCO+ could be destroyed 
during gas diffusion inwards, while most DCN survives. Such a scenario seems unlikely in 
light of recent model s of disk chemistry with diffusion that show HCN chemistry is faster 



than that for HCO"*" ( ISemenov fc Wiebell201ll ). A HCN deuterium fractionation mechanism 



that differs from DCO"*" chemistry is thus required to explain the observed DCN distribution. 
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Photochemistry can produce deuterium fractionation and iThi et al.l ( 120101 ) show that 
high HDO/H2O ratios in disk gas can be achieved through photochemistry at T > 100 K, 
where the HDO formation rate depends on the 0+HD reacti on. DCN can similarly form 
through CN+HD. This reactions has a barrier of ~3000 K (j Johnston fc BersohnI Il989l ). 
however, and formation of DCN from CN+HD seems unlikely to contribute significantly to 
the TW Hya DCN distribution on > 50 AU size scales in the absence of significant vibrational 
or electronic excitation of the CN. This is consistent with the HDO model results, where the 
same mechanism seems to be exclusive to the atmosphere of the inner disk where UV fluxes 
are very high. Since ultraviolet induced photochemistry should operate mainly in the disk 
atmosphere, resolved observations of multiple DCN transitions that can constrain the DCN 
excitation and vertical abundance distribution could be used to address the viability of this 
mechanism. 

Freeze-out onto grain surfaces has the potential to regulate the DCN distribution relative 
to DCO"*" because of the different binding energies of DCN and CO. Gas absorption onto 
grains is predicted to produce molecule-specific "snowlines" beyond which abundances drop 
dramatically. This mechanism clearly results in a centrally peaked abundance distribution, 
though warm deuterium chemistry still would be required to produce DCN in the first place. 
It seems unlikely, however, that freeze-out could regulate the HCN and DCN distributions 



on th e observed size scales. HCN has a comparable adsorption energy to H2O (l Aikawa et al. 
1996) . which, for the TW Hya disk, results in a snowline at a radius < 10 AU (lOi et al. 



20081 ) . This does not exclude the possibility that the DCN distribution is directly regulated 
by depletion onto grain surfaces at smaller radii. In fact, in the inner hot region of several 
disks, where HCN is detected at infrared wavelengths, the abundances are much higher than 
observed in outer disks at millimeter wavelen gths, suggestive o f an abundance drop at a few 
AU due to ice formation beyond this radius (jSalyk et al.ll201ll ). 



Beyond this snow line, ice forrnation products may be returned to the gas-phase through 



photodesorption (e.g. lOberg et al.ll2009[ ). Both laboratory experiments and observations of 



grain surface products show that deuterium fr actionation can be very efficient in interstellar 
ices (iNagaoka et al.l l2005l : iParise et al.l l2006l ). Theoretically the DCN/HCN ratio in ices 
can reach values of 2 x 10 ~^, similar to the ratios predicted from pure gas-phase chemistry 



( lAikawa fc HerbstI Il999l : IWillacyl l2007l ) . Recent determinations of low abundance s of cold 
H2O gas toward TW Hya and DM Tau (IBergin et al.ll2010l : iHogerheijde et al.ll201ll ) suggest, 
though, that H2O ice evaporation is not efficient outside of a few AU radius. Because of the 
low abundances of HCN and DCN relative to H2O, these species are likely to act as minor 
impurities when formed in the ice, and desorb only where H2O ice desorbs. Thus, while we 
cannot rule out ice photodesorption as a contributor to the observed DCN distribution, it is 
unlikely to be the dominant process. 
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Having dismissed dynamics, photochemistry, and mechanisms requiring grains, as plau- 
sible routes to significant deuterium fractionation in the outer disk, this leaves additional 
gas-phase formation routes of DCN as the main alter native to the r eactions involving CH2D+. 
In the gas-phase, DCN can form from HNC+D (e.g. IWillacyll2007l ). The deuterium fraction- 



ation in HCN depends then on the gas- phase D/H ratio. The D/H ratio is predicted to be 
< 10~^ at large disk radii (250 AU in a generic dis k model), whi l e it reaches and possibl y 
exceeds lO'^ in the inner disk midplane (17-50 AU) JwillacvlboOTi : Iwillacv fc Woodslboool ). 
Thus, reactions between D and HCN provide a plausible pathway to deuterium fractionation 
in HCN in some regions of the disk, and detailed modeling is needed to provide predictions 
on the vertical and radial abundance profiles resulting from this formation mechanism. 

In summary, several reaction and destruction mechanisms may contribute to the DCN 
abundance pattern observed in the TW Hya disk. However, warm CH2D"''-driven deuterium 
chemistry in the inner disk and cold H2D^-driven deuterium chemistry in the outer disk pro- 
vide are likely to dominate over the alternatives at the disk radii probed by the observations, 
and these pathways form consistent picture on their own. A warm CH2D"'"-driven deuterium 
chemistry was also the favored explanation for the observations of DCN in the Orion bar 
(jParise et al.ll2009l ). To determine definitively which formation pathway (s) dominate in disks 
will require interferometric observations with higher spatial resolution, observations of re- 
lated molecules and isotopologs, and more detailed models of DCN radial distributions for 
the different chemical scenarios. For example, if photodesorption is important, then DCN 
should spatially correlate with other grain surface products, and a powerful photochemistry 
may enhance the ^ ^N/^^N ratio in HCN, analogous to what is observed in Titan's atmosphere 
Jliang et al.ll2007l ). 



Whichever DCN formation pathway dominates, the substantial deuterium enhancement 
in warmer regions of the TW Hya disk challenges the conventional wisdom that the high 
deuterium fractions in corn ets is necessary evidence for a cold (T<30 K) chemical history 
(IMumma &: Charnleyll201ll ). Most comet deuterium fractions are measured from the D/H 
ratio in cometary water. In Oort cloud comets this ratio is ~ 3 x 10~^, enhanced by an order 
of magnitude compared to the cosmic D/H ratio. Higher D/H ratios of 0.0023—0.025 have 
been observed for DC N/HCN in Hale-Bopp, consistent with models of gas-phase deuteration 
followed by freeze-out (IMeier et al.lll998l : iBlake et al.lll999[ ). In light of the TW Hya results, 
even this high ratio does not imply a cold, T < 30 K origin of the comet material, but may 
instead reflect the existence a lukewarm pathway to deuterium enhancements in the Solar 
Nebula. The radial distributions of DCN and DCO"^ in disk material around TW Hya and 
other T Tauri stars should be revisited by ALMA at higher resolution to investigate whether 
the observed deuteration chemistry trends persist at smaller scales. Until then, care should 
be exercised in interpreting the D/H ratios in cometary material. 
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In general, resolved distributions of molecular emission in disks have the potential to 

put stronger and more direct constraints on key aspects of chemical evolution than (global) 
abundance ratios. In the case of DCN, the DCN/HCN abundance ratio of ~ 10~^ is consis- 
tent with a range of chemical models that assume different formation pathways. However, 
the observation of different radial distributions of DCN and DCO"'" immediately suggests the 
presence of different formation pathways for these species and rule out that DCN is mainly 
formed through reactions with H2D+. 
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Fig. 1. — Integrated images of the TW Hya disk in (left:) DCN J=3-2 e mission (ALMA 
science verification) and (right:) DCO"^ J=3-2 emission (SMA data from iQi et al.l l2008l ). 
The contour levels are at 50%, 75% and 90% of the peak value. The cross marks the peak 
of the continuum emission, which locates the position of the central star. 
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Fig. 2. — The best-fit radial column density distributions of DCN in the disk around TW 
Hya from 10 to 100 AU for a falling power law {a = —2, thick line), a rising power law 
(a = 2 dashed line) and a constant (a = dotted line). 



- 12 - 



4 

" 2 



[JCN 3-2 3 '"' 


3.51 


3.31 


3.11 


2.91 


2.71 


2.51 

, s 


. 2.31 

n r 


^ 2;.11 
★ 

n 


Model 1 3.71 

{a=-1) 

e 


3.51 


3.31 


3.11 


2.91 


2.71 


2.51 


2.31 


2.11 

★ 


Model 2 3.71 

(a=0) 

Ik- 


3.51 


3.31 


3.11 


2.91 


2.71 


2.51 


2.31 


2.11 

★ 


Model 3 3.71 

{a=2) 

★ 


3.51 




3.31 


3.11 


2.91 


2.71 


2.51 


2.31 

★ 


2.11 

★ 



4 2 0-2-4 
Aa(") 



Fig. 3. — DCN 3-2 velocity channel maps toward TW Hya. The top row shows the ob- 
servations, and subsequent rows show models corresponding to power law column density 
distributions with indices a = —1 (Model 1), (Model 2) and 2 (Model 3). The data is best 
matched by Model 1, as can be seen by comparing the shape of the central channels and the 
extent of extreme velocities. The contours are 0.01 x [2,4,6,8,10,14,18,22,26,30] Jy beam~^. 



